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We propose a new strategy to improve the self-assembly properties of DNA- 
functionalised colloids. The problem that we address is that DNA-functionalised 
colloids typically crystallize in a narrow temperature window, if at all. The underly- 
ing reason is the extreme sensitivity of DNA-mediated interactions to temperature 
or other physical control parameters. We propose to widen the window for colloidal 
crystallization by exploiting the competition between DNA linkages with different 
nucleotide sequences, which results in a temperature-dependent switching of the 
dominant bond type. Following such a strategy, we can decrease the temperature 
dependence of DNA-mediated self assembly to make systems that can crystallize in 
a wider temperature window than is possible with existing systems of DNA func- 
tionalised colloids. We report Monte Carlo simulations that show that the proposed 
strategy can indeed work in practice for real systems and specific, designable DNA 
sequences. Depending on the length ratio of the different DNA constructs, we find 
that the bond switching is either energetically driven (equal length or 'symmetric' 
DNA) or controlled by a combinatorial entropy gain ('asymmetric' DNA), which re- 
sults from the large number of possible binding partners for each DNA strand. We 
provide specific suggestions for the DNA sequences with which these effects can be 
achieved experimentally. 
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I. INTRODUCTION 



The high selectivity of hybridisation of DNA makes it an interesting molecule to be used 
as "smart glue" in the self-assembly of complex, nano-structured materials. Some of the 
advantages of DNA as a selective linker are that it can code for a large variety of specific 
linkages, it is readily available, it can be used under near ambient conditions and the bind- 
ing process is reversible. Consequently, applications of DNA-mediated self-assembly range 
rom computational biology [1] and the assembly of scaffolded "DNA-origami" structures 
^ to the development of targeting strategies (see e.g. 3-5]). Much experimental work has 



focused on the application o 
complex colloidal structures 



' DNA as a selective linker that enables the self-assembly of 
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25]. Several examples have been reported in the literature 



of (relatively simple) DNA-linked colloidal crystals, either consisting of nano colloids (see 



e.g 
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16]) or of micron-size colloids with complementary DNA coating (see e.g. [11]). 
In parallel, theoretical investigations have provided insight into the factors that influence 



DNA-mediated self-assembly |26l-l30|. 

There are, however, some disadvantages associated with the use of DNA as a tool to link 
(nano)colloids. In particular, the very factors that lead to the specificity and reversibility 
of DNA linkages also cause the strength of DNA mediated interactions to depend strongly 
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-|23| • This 



on the external conditions, such as the temperature or ionic strength 
sensitivity results in an abrupt onset of aggregation as the temperature is lowered - a phe- 
nomenon that can be problematic for self-assembly because it narrows the "window" of 
conditions within which reversible self-assembly of ordered structures is possible. If this 
window is missed, self-assembly either does not take place at all or results only in disor- 
dered aggregates. It is therefore important to explore possible approaches to combine the 
high selectivity of DNA-mediated interactions with a more gradual response to external 
conditions. 

In the present paper, we use numerical simulations to demonstrate that 'competing' 
DNA interactions can be used to create colloidal systems with a more gradual temperature 
response. Specifically, we consider a binary mixture of colloids (X and X') that, unlike most 
systems studied so far, are functionalised by a mixture of different DNA strands (a, (3 and 
a?, (3', respectively). These DNA sequences are chosen such that a can bind to both a' and 
f} 1 (and similarly a 1 to a and /?), but cannot bind to j3' (we will show that such DNA 
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sequences can be readily designed). Importantly, the a-fi' and a'-(3 linkages are weaker 
than those between a and a'. As a result, a-a' can form at a higher temperature than 
the linkages involving or /3'. However, whereas a given a, a' pair can form only one Di- 
et' linkage, it can form two weaker linkages (ct—fl' and a'-/3). This leads to a competition 
between the different types of linkages in which both their free energy of formation and 
combinatorial entropy effects,which depend on the number of different ways in which the 
inkages can form, play a role - as is characteristic for systems with multivalent binding 
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37). Here, we investigate how the majority of the linkages can switch from one type to 



the other, not only as a function of the temperature and the difference in binding strength 
of the two linkage types, but also as a function of the surface coverage and length ratio of 
the DNA constructs. 

The remainder of this paper is organised as follows: in Sec. [TT] we introduce the model 
and the Monte Carlo algorithm. Sec. IIHI reports our results. In Sec. IIII Al we consider the 
symmetric model, in which the length of the competing DNA constructs is the same (L a = 
Lp). We show how the effective pair interaction depends on the hybridisation free energy of 
the DNA sticky ends and how the energetically driven switching from one linkage type to 
another broadens the association-dissociation transition of the particles; we rationalise these 
findings in the context of a Mean-Field model. In Sec. IIIIBI we consider the asymmetric 
model (L a < Lp). We show how the length ratio of the different DNA constructs can be 
used to enhance the bond switching through combinatorial entropy effects. Finally, in Sec. 
HV] we indicate how our approach could be implemented experimentally. 



II. THE MODEL 



For the sake of simplicity, we use the relatively simple but well-tested model of ref. [31 1 
to describe the interaction between DNA-functionalised colloids. 

Following ref. 31], colloids are coated with double-stranded (ds) DNA that is shorter 



than its persistence length of i 50 nm 38) (corresponding to 150 nucleotides 39) ). The 



dsDNA is terminated by a short single-stranded (ss) DNA sequence ('sticky end') that can 
hybridise to complementary ssDNA. Binding is only allowed between colloids of type X and 
X' . DNA is modelled as thin rods (see Fig. [1]) randomly tethered to the colloidal surface. 
For typical surface coverages, the steric repulsion due to self-avoidance of the DNA strands 
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FIG. 1: Model system investigated in the present work. Pairs of colloids (X and X') are func- 
tionalised by two families of short DNA fragments (a, f3 and a', ft'). Colloids feature two kinds 
of linkages: a-a' (a) and a—/3'+a'—/3 (b). When hybridised the conformational space of rods re- 
duces to QaP'/aa'- Due to the excluded volume interaction between colloids and dsDNA, available 
configurations are also reduced when two colloids approach each other (c and d). 

(2 nm diameter) is negligible compared to the entropic effects involved in the binding of 



tethered strands [3l|] and we therefore treat the DNA as non-self-avoiding. Except for the 
DNA funtionalization, the colloids are assumed to be smooth and hard. We consider two 
types of DNA (a, (3 on X and a', (3' on X', Fig. [I]) denned by different ssDNA end sequences, 
and (later on) different lengths L a , Lp. The characteristic distance between DNA strands on 
a given colloid is S = a/ 'A/(N a + Np), where A is the total surface area of that colloid (say 
type X) and N a and Np are the total number of a and ft strands on X (similar expressions 
apply to X'). 
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Here we focus on the interaction between two parallel, planar surfaces. The pair potential 
between spherical colloids (of radius R) can be computed from these planar surface interac- 
tions via the Derjaguin approximation 40] . In Ref. 31] we have shown that this is a reliable 
approximation when R/L > 10. From geometry, it further follows that the interactions are 
strictly pairwise additive if R/L > 6.5, assuming that the hard cores of the particles can 
come into contact. For smaller particles, curvature and so-called three-body effects may 
cause deviations of the exact interactions from what is predicted here. Nevertheless, we 
expect that bond switching can also occur in these systems, albeit at different values for the 
relevant parameters. 

In the simulations, we consider two square planes with side L = \J~A and periodic bound- 
ary conditions in the directions parallel to the planes. We fix L a + Lp = 40 nm [31] . In 
real units, the width of the simulation box is L w 0.5 fim. Varying N a and Np, we obtain 
characteristic inter-chain separations S = 0.53, 0.75, 1.06 • L a . These values span typical ex- 



perimental values, e.g. 10j,llll]- We verified that averaging over different random realisations 



of tethering points does not alter our results within a 2% of tolerance. 

The hybridisation of pairs of ssDNA depends on the specific nucleotide sequences of the 
individual strands and on the solvent properties (temperature and salt concentration). For 
the present paper it is important that ssDNA sequences can be designed in such a way that 
only the following linkages are possible: a-a', but also a-f3' and a'-(3 (see Fig. [T]a and b). 
In Sec. IIVI we show that it is possible to design nucleotide sequences that will yield this 
behaviour. 

We denote the hybridisation free energy of two free ssDNA by AGq x , where x = a 
for an ot-ot' linkage and x = (3 or x = (3' for an a—/3' or a'— ft linkage. We assume that 
AG 0a < AG p> = AG p. The binding free energy between two surface-bound DNA strands 



must also include a configurational entropy term [31], |37|] , which accounts for the reduced 
freedom of motion of the strands upon binding (Fig. [1]). If the dsDNA can swivel freely on 
the surface and the ssDNA is connected flexibly to the dsDNA [l^] , then the configurational 
space of the DNA rods is a circle (Q aa > and Q a p' in Fig. CD), that is bounded by the two 
surfaces. In App. IA1 we give an explicit expression for Q a p> as a function of L a and Lp. As 
can be seen from Figs. [T]c and d, fl a and Qp are equal to 2nL 2 a jp if the distance between the 
planes h is larger than L a /p and equal to 27rL a /ph if h < L a /p. As a result (see 0] and 
App. [B]) the hybridisation free energy of two tethered strands AG ax (x = a' or (3') is given 
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by 

AG ax = AG Qx -k B T\n(--^- ) , (1) 

where po is the number density corresponding to one molar p Q = 6.022- 10 23 /liter. r is a non- 
universal factor that depends on the details of the linkage formed between the two dsDNA 
strands. As explained in App. [Bj slightly different expressions for r result depending on 
how the sticky ends are coarse grained. However, the resulting differences in the predicted 
binding strengths are tiny and therefore irrelevant for the present discussion. 



In each MC move 



311 ]. we randomly select one of the atj (or a',) functional arms and 



attempt to make, break or switch a linkage. First we consider the list of all the free x'-, 
(x = a, P) that can bind to atj, possibly including the partner to which otj is already bound. 
A linkage aj- x' q , is created with probability p q > = exp[— 0AG ajX > r ]/Qj. The probability that 
no linkage is formed is given by po = l/Qj , where Qj = 1 + J2 X ' ex P[~ (3AG a v,]. AG ax ' is 
computed using Eq. ([T|). In order to enhance the sampling of the model we also implemented 
a 'linkage-swapping' MC move that attempts to switch between a single a-a' linkage and 
two weaker a-(3'+a'—fl linkages. Details are given in App. O 



III. RESULTS 



A. Symmetric DNA model (L a = Lg) 

We first consider systems in which the competing DNA constructs have equal length L a = 
Lp. Fig. [2]shows n a , the average number of a-a' linkages, and np, which denotes the number 
of a-/3'+a'-/3 linkages. n a and np are plotted as a function of AG 0q . There is a linear relation 
between the hybridisation free energy and the temperature AG 0a = AH Qa — TAS 0a , where 
the hybridisation enthalpy/entropy (AH 0a / AS 0a ) are negative constants. Hence decreasing 
AGo a corresponds to a decrease in T. For real DNA sequences, 5 AG = AGop — AGo Q also 
depends on T (see Sec. HV|) . However for the sake of simplicity in Fig. [2] we sketch our results 
at constant SAG. 

Upon decreasing AGo a from an initial situation where all DNAs are unbound, the stronger 
a-a' linkages form first. For the specific choice of simulation parameters listed in Fig. &\ this 
happens when AG 0a m —5k B T. As we lower AG 0a even more, the n a linkages disappear in 
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FIG. 2: Number of a— a! linkages (black curves) and of a—ft'+a'—ft linkages (coloured curves) 
versus AGo a (taken in unit of ksT). In each simulation 5 AG = AGop — AGo a is kept constant 
and equal to the value reported at the crossing point (n a = np). Full lines refer to MC simulations 
that only use single bond rearrangements, while dotted lines refer to simulations that also employ 
'linkage-swapping' moves (App.[C]). We have used N a = Np and S = 0.75L a , while the distance 
between colloids is h = L a = Lp. 

favour of the weaker bonds between a— ft' and a'- ft. This may seem counter-intuitive at first, 
but the reason is simple: replacing a single a-a! linkage by two weaker ones (a-ft'+a'-ft) is 
favorable as long as the total gain in binding free energy upon forming these two linkages 
outweighs the loss in binding free energy when breaking an a— a' linkage. The switching 
between the different bond types in Fig. [2] is driven by the difference in the hybridization free 
energies of the individual linkages and the main effect of increasing 5 AG is to decrease the 
AGo a value of the crossing point where np = n a . As will be shown below, the combinatorial 



entropy of the overall system 3l|, |37| plays a less important role in this particular case. 

The behaviour shown in Fig. [2] can be understood in a more quantitative way using a 
mean-field theory. We define s a (h) and sp(h) as the average configurational entropy cost of 
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FIG. 3: Monte Carlo simulations (full lines), Mean Field Eq. [4] (broken lines), and self-consistent 
Mean Field Eq. [5] results (dotted lines) for n a and np. Part (a) reports predictions for the single 
linkage model with S a = 1.06 (first column), S a = 0.75 (second column), S a = 0.53 (third column), 
h = L a (first row), and h = 1.5L a (second row). Part (b) reports predictions for the competing 
linkages model with simulation parameters like in Fig. [2j 



an ol-ol' and at—/3' linkage 



s a (h) = — I dx' exp 



1 a J A, 



AS^(L a ,L c 



sp{h) = -7-1 dx' exp 



AS^(L a ,L fi )i 
k B T 



(2) 



where A a and Ap are the circles on X' enclosing all the possible a' and /3' (tethered at position 
x') which could be hybridised by a on X-DNA coated colloids represent a multivalent 
binding system in which each DNA strand can form several different linkages. In the low 
binding regime (n a + ng/2 N a ) this combinatorial entropy contribution is accounted for 
by noticing that the number of possible a— a' and a-/3' linkages are A a S~ 2 and ApS^ 2 , where 
S a {Sp) is the mean distance between a ((3) strands. It follows that the partition function 
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Z and the average fraction of the two linkages (p a = n a /N pp = np/N) is: 

'N a \ f2N a - 2n a \ /s a (h)e-P AGo «\nc, fsp{h)e- pAG ^\^ 



2n a +np<2N a 



Up J\ SI/A, J V S 2 p/Ap 



l ^ s a {h)e^ AG ^ + 2 sp(h)e~ 0AG w + fSp(h)e^ AG ^\^ 



Pa(S a , Sp) 



Sl/A a ~ Sj/Ap V Sj/Ap 

s a {h)e-P AGo « 
ZV*«Sl/A a 

2sp(h)e-P AG ^ +2(sp(h)e-P AG A 2 /(S 2 p/Ap) 



(3) 



Pp(S a ,Sp) - zV N «S 2 p/Ap ' (4) 

Fig. [3] compares the mean- field predictions based on Eq. (Tj0) for the single linkage model (mI 
(part a) and for the competing linkages model of Sec. [TTl (part b) with the results of the MC 
simulations. There is qualitative, though not quantitative agreement between theory and 
simulations. 

A better agreement is obtained when we improve our estimate of the combinatorial en- 
tropy terms. If n a linkages are already present, the number of ways an extra linkage can be 
added is smaller than if n a = 0. This is, because the mean distance between un-hybridised 



a' strands increases like S a / y 1 — n a /N (and similarly for Sp). Instead of exactly dealing 
with this correction we used Eq. (JJ]), while correcting S a and Sp by the average number of 
linkages (n a and np) which are then computed self-consistently in the following way 
n a _ / S a Sp \ n£ _ / S a Sp \ 

(5) 

Dotted curves of Fig. |3] show the predictions based on the solution of Eq. [5j As can be 
seen from this figure, the agreement between theory and simulation is now satisfactory. 
Although producing different profiles, Fig. |3b shows that the two different estimates of the 
combinatorial entropy (Eqs. H]and|5]), place the bond switching transition (where n a = np) 
at the same AG 0a . Indeed, using Eq. (j4j) without any combinatorial prefactors (but allowing 
for each a /a' strand only one a! /a and one partner) we find that AGoa at the bond 
switching transition decreases by J$ 3ksT. This indicates that the bond switching observed 
here is mainly an 'energetic' effect, having to do with the different hybridization free energies 
of the individual bonds. 

The fact that the dominant bond-type between colloids switches as the temperature is 
decreased, gives rise to a more gradual temperature dependence of the attractive interaction 
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FIG. 4: Minima of the pair potential per unit area as a function of (a) AG$p and (b) the fraction 
of hybridised linkages. We report results for systems with only np or n a linkages (AGo a = 00 
or AG0/3 = 00), and for the competing linkages model with three values of 5 AG. We have used 
S = 0.75 L a , L a = Lp = 20 nm, and N a /Np = 3/7. 



between the colloids. To see this, we consider the effective pair potential F(h), where h is 
the separation between the two surfaces. Following Ref. 31] ; F(h) can be computed using 
thermodynamic integration starting from the case where the DNAs do not hybridise. In 
that case the pair potential is purely repulsive and is given per unit area by: 



Frep(h) = T^-max 



+ — 2 max 



(6) 



If we now consider colloids coated wit 
can generalise the arguments of Ref. 



1 two different types of linkers with constant SAG, we 



F(h) 



Fe P (h) + - 



3jj to show that the full potential is given by 

/•AGoa 

/ dAGd Q < n a + np >AG' 0a ,6G ■ 
J 00 



(7) 



S a and Sp are the mean spacing of the a and f3 strands (S a = a/ 'A/N a , Sp = v/ A/Np) 
and A denotes the unit area inside which linkages are counted. In Eq. (J7|) the average is 
taken with AG$p = AGq q + 5 AG. For the lower limit of the integration, we choose a value 
for AG' Qa such that < n a + rip 1. Fig. H]a shows the dependence of the depth of the 
attractive well of F(h) (usually located near h = L due to the steric DNA-surface repulsion 



3jJ ) as a function of AGop for three different values of 5 AG. For sake of comparison we also 
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show the results for a system with the same distribution of a and (3 sticky ends, but in which 
the a-a! linkages are forbidden (equivalent to taking AGo a = oo). This choice allows us to 
compare between a competing and a single linkage model with equal maximum number of 
possible bonds. At low AG0/3 only np linkages are present, and all the curves approach the 
same binding strength. Importantly, the figure illustrates that the strength of the attraction 
changes more gradually with temperature when competing DNA linkers are present - this is 
especially evident at high SAG. Experimentally, the more gradual response to temperature 
means that the DNA coated colloids display a broader association-dissociation transition and 
have a larger range of conditions under which they can form ordered assemblies (instead of 
kinetically disordered aggregates). 

Fig. H|> shows F min versus the fraction of hybridised bonds / = (n a + np)/N tQt , where iVtot 
is the maximum number of available linkages. Compared to single linkage models (black 
curves), the competing linkages model acquires a reasonable attractive well at lower /. It is 
worth remembering that the strong, high temperature a— a! linkages are replaced with weaker 
a-/3' linkages as the temperature is lowered. The fact that the a-a' linkages disappear before 
they become prohibitively strong means that kinetically trapped configurations should be 
automatically avoided, which represents a great experimental benefit. Indeed the formation 
of an 'irreversible' a— a' linkage starting from two 'dynamic' a—fl'+a'—fl ones needs the break- 
ing of two independent linkages. Under bond-switching conditions, this requires an average 
time (At ~ exp[— 2(3AG a p>]) even longer than the life-time of at— a! (At ~ exp[— 0AG aa i]). 
Nevertheless, Fig. El does show signs of equilibration problems. Using only sequential single 
linkage MC moves (full lines in Fig. E} we found that, at low AGo a and high 5 AG, the 
system was not able to completely remove the a— a' linkages. The problem is due to the 
fact that although the formation of two weaker bonds is thermodynamically favourable the 
system needs to first break a strong linkage. To equilibrate the system (dotted lines in Fig. 
EJ) it was necessary to use the MC move described in App. [Cj which swaps between a strong 
linkage and two favourable weaker ones. This implies that, in order to engineer kinetically 
accessible experiments, it is important to design nucleotide sequences with a small 5 AG, and 
to think about possible strategies to push the bond switching transition to higher AGo a . 
This will be investigated in the next section. 
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B. Asymmetric DNA model (L a < Lp) 

In the previous section we have discussed symmetric systems (L a — Lp) for which the 
switching from one strong bond to two weak bonds was mainly driven by the gain in binding 
free energy, due to the difference in hybridization free energy of the different linkage types. 
However, the bond switching can be further enhanced and, more importantly, kinetically 
enabled, by modifying the number of possible binding partners that each type of linker 'sees' 
on the opposing surface. In doing so, we take full advantage of the multivalent nature of the 
system, going from a mostly 'energetically' driven bond switching mechanism to one that is 
driven by combinatorial entropy effects. 

Eqs. Hlshow that multivalence controls the appearance of a-a' or a—P'+a'-ft linkages by 
A a / S 2 and Ap/Sp, where A a and Ap is the area of the surface that enclose the tethering 
points of all the possible targets of a (a' and /?'), while S a /p is the mean distance between 
DNA strands. Using simple algebra we find 



Aa_ = ^L 2 a N a 
S 2 A 



h 2 



Ap _ n(L a + LpfNp 



h 2 



(2L Q ) 2 J S 2 A L (L Q + Lp) 2 . 

where we have introduced the number of a//3 strands (N a /p) for unit area A. One obvious 
way to favour the weaker a-/3' over the stronger a-a' bonds is to decrease the concentra- 
tion of a-terminated DNA (N a ) relative to that of the more weakly binding /3-DNA (Np). 
However Eq. (JSJ) suggests that a-a! linkages may be suppressed even more efficiently by 
decreasing the length L a . For instance, at h — L a , the same effect of having a relative 
concentration of a strands equal to 1%, can be obtained by an asymmetric system with 
L a = l/3Lp. Fig. |5]a shows n a and np under the same conditions as those considered in 
Fig. |2j but using L a j ' Lp = 3/17 (L a + Lp = 40 nm). As compared to Fig. [2] the AG 0a value 
for which n a = np is ~ 3ksT higher if 5AG = 3ksT, m 8ksT higher if 5AG = 8ksT, and 
~ 14/c^T higher if 5AG = lAksT. We used both local and linkage swapping MC moves 
(App. [Ujl . Interestingly, the n a linkages first saturate at a fraction well below 0.5 (= N Q /N). 
This means that many a and a' sticky ends remain free and that the np linkages are formed 
much earlier than in the symmetric case (see Fig. [2]). The plateau value where n a saturates 
depends on the DNA mean distance S (Eq. [8]) and ranges from 0.14 if S = 1.06 to 0.31 if 
S = 0.53. The value of AG 0a where n a = np also has a dependence on the coverage density, 
when using two lengths for the DNA constructs, but this dependence is found to be small 
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FIG. 5: (a) Same as in Fig. Hi (N a = Np) except that L Q = 0.3, Lp = 1.7 and h = 0.3. (b) 
Colloidal pair potentials F(h) (Eq. EJ) at six different AGo a and 5 AG = 8ksT. The vertical line is 
the distance at which plots of part (a) have been recorded. 

(;$ 2ksT for the tested cases). 

Fig. |5]& shows the effective pair potential per unit area at different values of AG a, keep- 
ing AGop — AGo a fixed at 8k bT. For large values of AGo a , there is a local minimum in 
the potential at h = L a . This minimum becomes a global minimum as AGo a is lowered. 



Important features of Fig. |5]fr are the presence of repu 



low AGqo. By using a Derjaguin approximation 
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sive tails that become shoulders at 



4fJ it is possible to design sensible 



interactions that are useful to stabilise complex assemblies. 



IV. DESIGNING THE STICKY END SEQUENCES 

The analysis described in the previous sections allows us to predict the relation between 
the formation of a— a', a— ft' or a'— ft bonds and the hybridisation free energies of the indi- 
vidual strands. In addition, we can now understand the effect of surface coverage and chain 
length. We find that the dsDNA lengths (L Q and Lp) are the most relevant parameters that 
control the bond switching transition. 

We now consider possible choices for the sticky end sequences that in experiments would 
result in a broadened dissociation transition of the colloids. We need to design our ssDNA 
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FIG. 6: Two possible structures of the reactive end fragments (thick lines) with hybridised states 
that involve different subdomains (coloured segments labeled by a, b, and c). a* , b* and c* are the 
complementary strands. 

sequences such that the a sequence binds more strongly to a' than to /3' (and similarly 
for a'). If we limit the discussion to Watson-Crick base pairs, then different domains of a 
hybridise with a' and with fl'. In Fig. [6] we consider two architectures for the ssDNA in 
which the a-a! linkage involves (a) the full length of the reactive fragments or (b) pieces of 
ssDNA which are not used in the weak bonds. As a proof of concept in Tab. HJ we report 
possible sequences for these two families (A for Fig. |6]a and B for Fig. [6f>). 

Sequences of type A have been chosen with a = c (Fig. [6]a). This straightforwardly 



balances the hybridisation free energy of the ct—/3' and the a'-/3 linkages 42( . The bottleneck 
is the possibility of fi—fi' linkages, weaker than a—/3'. This problem is avoided by using the 
scheme of Fig. [6]& and less degenerate sequences {B\ and B2) to further enforce the selectivity 
of each subdomain. In particular B2 has been assembled using for the weak linkages two 
couples of sequences with nearly equal hybridisation free energy |4l| . 

As input for our model we use the hybridisation free energies of ssDNA fragments in 



solution. We used the "DINAMelt' estimates of the ssDNA binding free energies 43|. For 
more details on the procedure, see Appendix [B] Fig. [7] shows the relation between AGo a and 
AGqq — AG0/3 that is predicted for ^(n) and ^(n) listed in Tab. Hand different values of n. 
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seq. Ai(n) 


a = 5' - [C n ][T n ][C n ] - 3' 
P = &-[C n ][T n ]-& 


a' = 5'-[G n ][A n ][G n ]-3> 
P' = 5' - [G n ][A n ] - 3' 


seq. A2{n) 


a = 5'-[r n ][c n ][r n ]-3' 

P = 5'- [T n ][C n }-3> 


a' = 5' - [A n ][G n ][A n ] - 3' 
P' = 5' -[A n ][G n }-3' 


seq. B\(n) 


a = 5' — [(TCT) n ][(GTG)n] — 3' 
/3 = 5' - [{TGT) n } - 3' 


a' = 5'- [(ACA) n ][(CAC) n ] -3' 
P' = 5' - [(AGA) n ] - 3' 


seq. B 2 (n) 


a = 5' - [TTGAGAAATCC] [C n ] - 3' 
P = &- [AAGATTGATCC] - 3' 


a' = 5' - [G n ][GGATCAArCrT] - 3' 
/3' = 5' - [GGATTTCTCAA] - 3' 



TABLE I: Nucleotide sequences giving rise to the hybridised states reported in Figs. 0a and [6]6 
(i?2 partially taken from [4l|). Brackets group different subdomains (a, b and c) as defined in Fig. 
[6j while subscripts stand for repeated nucleotides. 



Based on the MC results of the previous section, we can predict the different bonding regimes 
both for DNA strands of equal length (Fig. 0a) and for the asymmetric case L a f 'Lp = 3/17 
(Pig.©). 

Fig. [7] can be used as the starting point to design possible experiments. First, given 
a certain colloidal architecture, MC results allow to draw the region where none, n a or rip 
linkages are expected. For given ssDNA sequences, it is then possible to predict the transition 
temperatures computing AG 0a and SAG and overlapping them with MC results like those 
in Fig. [71 Sequences can be optimised (e.g. changing n in Tab. [I]) to avoid kinetically trapped 
configurations. 

V. CONCLUSIONS 

In this paper we have proposed and tested a strategy that could make the crystallization 
transition of DNA-functionalised colloids less sensitive to external conditions. 

Specifically, we have considered a binary system of colloids (X and X') covered by two 
families of reactive sticky ends (a, P on X and a', P' on X'). Exploiting the selectivity 
of DNA, we have shown that it is possible to design sequences that only allow for binding 
between a-a', a~P' and a'-p. We choose the hybridisation free energy of a and a' in 
solutions stronger than the other two possible pairings with equal hybridisation free energy. 
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FIG. 7: Full black lines are MC results which outline the kind of linkages (none, a-a' , or a-(3') 
in the {AGo a , AGo a — AG0/3} plane. Full circles and open squares differ in the way the reactive 
sticky ends are coarse grained (see App.[B|). Part (a) refers to the symmetric system {L a = Lp, 
h = L a ), while part (b) to L a = 3/17Lp, h = L a . Broken lines list some sequences of Tab. [J 
within a temperature range (recorded every 10 °C) where the crossover between different linkages 
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Because a or a' participate in all the possible linkages, the strong a-a' linkages compete 
with the weaker ones (a-/3'+a'-j3). However, the number of possible weaker linkages is twice 
the possible number of the strong ones. 

We have demonstrated that a bond switching transition is possible: upon decreasing 
the temperature the first linkages to appear are the strongest a-a' ones while at lower 
temperatures the system can gain free energy by replacing a-a' with a-f3'+a'-(3. For 
symmetric systems (in which the length of the DNA constructs is equal), the bond switching 
is energetically driven. This means that the hybridisation free energy of having two weaker 
linkages is lower than the hybridisation free energy of a single strong one. When choosing 
the length of the a and a' linkers shorter than /3 and /3' (asymmetric model) the bond 
switching transition is enhanced and occurs at higher temperatures than in the symmetric 
case. Here the transition is driven by the combinatorial entropy gain related to the fact that 
an a strand, for instance, can bind more j3' than a'. 

The main effect of the competing DNA linkages is that the resulting effective inter-colloid 
pair potential is less strongly temperature dependent than is observed in the conventional 
case where only a single type of linkage is possible. This enhances the experimental control 
over the self-assembly of DNA functionalized colloids. A further advantage of the proposed 
strategy is that the strong linkages are replaced as temperature is lowered, which could be 
used in in step-wise assembly schemes. 

Our procedure predicts the temperature range where transitions between different kinds of 
linkages are to be expected in experiments. Moreover, we have shown how one can optimise 
the DNA sequences and colloid architectures such that the linkage switching transition 
remains kinetically accessible. 
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FIG. 8: Two possible ways by which the colloids (broken lines) can cut the configurational space 
of the hybridised rods £l a p. Note that we choose a reference frame with the distance between the 
surface attachment points of the sticky ends d taken along the vertical direction. 

Appendix A: Computation of £l a p 

Here, we derive expressions for the configurational space of the hybridised strands, Q a p, 
(Eq. [1]) as a function of the distance between the tethering points of the hybridised sticky 
ends d, the rods lengths L a , Lp (L a < Lp), and the angle 9 between d and the vector normal 
to the colloid surface (see Fig. E}. If 7« is the angle 1 between L a and d, depending on the 
colloidal distance and the rod lengths, it may be 7 a < tt/2 or 7 Q > n/2. In the first case 
(Fig. Eb) & a p may have none (Eq. lAlfl . one (Eq. IA2I) or two cuts (Eq. IA3I) 

71 

Vt a p = 2nL a sin 7 a if 9 < - - j a (Al) 

7T 7T 

tt a p = [2tt - A(/)(7 a ,6»)]L a sin7 a if - - 7 Q < 9 < - - jp (A2) 

7T 

Q a p = [2tt - A<p(-f a ,9) - Aif(jp,6)]L a sin-f a if - - jp < 9 (A3) 

where we have defined A<£>(7, 9) as the planar angle of a cone (of amplitude 7) which is cut 
by a plane tilted by an angle 9 

Aip(j,9) = 2ArcCos ! 



tan 7 tan 9 



1 "f a is given by cos7 Q = (d 2 + L' a — 2dLp)/2. A similar relation holds for -frj. 
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If 7 a > 7r/2 (Fig. [5]&) there might be a single cut by the plane where L a is tethered (Eq. 
IA4I) . or the second plane could also cut f2 a/ 3 (Eq. IA5I) 

7T 7T 

fi a/3 = Aip(n - 7 Q , 6>)L Q sin7 a if 9 > j a - - and 9 < - - 73 (A4) 

71 71 

Q a p = [Ay?(vr - 7 a ,0) - A(y9(7 / 3,6')]L a sin7 a if > 7 a - - and 9 > - - 73 . 

(A5) 

Notice that Q a p = when 9 < 7 a — 7r/2. This happens only if L a < Lp, implying that for 
the hybridization of asymmetric linkages 9 needs to stay inside a narrower region than when 
L a — Lip. 



Appendix B: Estimate of the configurational entropy 

Following 



31], we first consider ssDNA fragments in solution (Y, Z and YZ) in thermo- 
dynamic equilibrium Y + Z f± YZ. In the ideal limit, the partition functions of the three 
species are 

^y, 2^ and are the contributions to the partition function due to the internal degrees 
of freedom [31], A Y /z are the de-Broglie wave lengths, V and N y /z/yz the volume and 
the number of molecules, while v is related to the specific potential (Vb on d) by which the 
two fragments are hybridised. Using classic pair potentials would give, for instance, vq = 
(27r/c£T//c#) 3//2 in the case of a harmonic potential with spring constant kjj, and Vq = A/Sttw 3 
for a square-like potential with amplitude w. At equilibrium the relation hyz = Hy+Hz links 
the ratio of the internal partition functions (zy/z/Yz) with the number densities (py/z/yz) 
which can be expressed in terms of the equilibrium constant K = exp(— f3AG Y z) 

ZYZVQ = Pyz_ = }_ e -^G° YZ (B2) 
zyZz PyPz Po 

We now consider the case of tethered rods. When un-hybridised the partition func- 
tions of the free fragments are Qy = flyZy/A Y and Qz = VL z zz/A 2 Z) where fly/z is the 
configurational space available to the sticky ends (Fig. [p. When hybridised 

Qyz = \2\2 / dx i / dx 2 exp[-/3 Vb ond (xi -x 2 )] , (B3) 
a y a z Jtoi Jn 2 
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where Xi and x 2 are the coordinates of the two sticky ends. Taking the large kn limit (for 
spring potentials) or the small w limit (for square-well potentials) we find 

A y A| sm(7y + 7 Z ) 

where and 7^ are the angles between the rods and the vector linking the tethering points 
(d in Fig. E]), and flyz as defined in Fig. CD Finally using Eqs. (1B4j) and ( 1B2I) we can compute 
AG YZ 

AG Y z = In ( = -AG YZ + In ( L -i . (B5) 

It is important to note that the 1/ sin(7y + 7^) term which appears in Eq. (IB5I) is specific to 
the case in which the sticky ends are modelled as point particles, but that it has no physical 
meaning. For the purposes of the present work, different prefactors give tiny differences 
which do not affect any of the presented results (compare black symbols in Fig. [7]). 

The hybridisation free energies of the sticky ends in solution AG° XY (AG 0a and AG 0/ g in 
the main text) have been computed as reported in 43[ and implemented on the DINAMelt 
server 44j using salt concentrations equal to [iVa + ] = 60 mM and [Mg ++ ] = mM. 



Appendix C: Linkage swapping MC moves 

Here, we discuss a Monte Carlo scheme which switches between a strong linkage a— a' 
and two weak ones a-(3'+a'-f3 (see Fig. [9]). In the en-en' — > a-[3'+a'-(3 move, we randomly 
select an a 3 — a'j, linkage and two unhybridised (3 strands ((3' ql , (3 q ) from the B(j) B'(j') that 
could be connected to aj and a'j, (Fig. [9]). For the reverse move a-(3'+a'-(3 — >■ a-a! we use 
one of two similar schemes: A) We randomly choose two hybridised couples otj-f3' q , + a!,,-j5 q 
within the possible A coup i c with atj and a'y neighbours or B) we randomly choose an aj-(3' ql 
(a'j,-f3q) hybridised pair from the n^i (ri/3',2) set (np = np^ + np/p) and, subsequently, an 
a f-/3q ( a j-(3' q i) from the L(j) (L'(f)) set with acj and a'-/ that could be connected. In the 
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FIG. 9: (a) Implementation of a MC move between competing linkages: a— a! (left) and a-f3'+a'-/3 
(right), (b) B(j) is defined to be the number of unhybridised (3' strands that could be reached by 
ay. L(j) is the number of a'-, which are hybridised by a (3 and that could be reached by ay. Similar 
definitions hold for L'(j') and B'(j'). 



two cases the acceptance rules are given by: 

n a B(j)B'(j')e 



accA(aja'j' — > atjfl'qt + a'j>/3 q ) = min[l 



decuple + L{j) + L'(j') + 1 

acc A {a 3 (3' q , + a' f P q -> «,«;,) = min[l, ^1 + + + 1] 1 ( C1 ) 

S(j)S'(/)e~^ [AG(a ^8' )+AG! K'^ ) - AG ^' a i'^ 



accs(a J -a'v — > otj/3' q , + a4/3 g ) = min[l, — 



acc B (aj/3' ql + a' jt ^ q ->■ ayay-/) = min[l, — — — lUD ^ ; 1irp ,^ A | n ] (C2) 



' ' ' " ' " L K 1 + l)[L(j) + l] 

n«iL(j> flAG(ai ^ )+AG(a i' ft) " AG(ajt "i' )1 
' (n a + l)[B(j) + l][B'(f) + l} 

In algorithm I? we randomly decide to extract first an a or a' hybridised strand. 

Alternatively, we can bias the algorithm to choose couples with low configurational en- 
tropy cost. In the a-a' —> a-(3'+a'-(3 move, we randomly choose an a-a' linkage, while j3', 
is selected with probability p q i 

P«* = w ( , " W B (j)= Yl ex P [-PAS^\ aj P' g ,)] . (C3) 

B{3) P' q ,eBU) 
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Similarly, (3 q is taken with probability p q , where p q and W B (f) are defined as above flC3[) . 
In the reverse move, like in algorithm B, we first randomly select an a.j^' q , (a'j,(3 q ) linkage 
and then one between aZj,(3 q ( ctj/3' q i) with probability py (pj), and: 

exp[-/3AS (cnf) («,•«',)] ^ , n 

Pi' = ^ L(j) J Wl(j)= £ exp[-/3AS( cnf )K-4)] (C4) 

"j/ei(i) 

(similar definitions follow for pj and W' L (j')). Acceptance rules are then given by 

n a W B (j)W B (f)e~^ 2AG w- AGo «\ 



acc bias (aja' j i ->■ <x,/^, + q;^/3 9 ) = min[l, 



(n p>1 + l)W L (j) 



acc hias (ccjl3 a , + av/3g ->■ a,-a,-,) = min 1, — ! — ^ ] , (C5) 

" 3 q 3 3 (n a + l)W B (j)W B (j>) 

where we have defined 

WbU) = WbU) + exp[-/3A5( cnf )(a,/3;,)] W L (j) = W L (j) + exp[-(3AS^ { \a 3 a ! f )} 

(C6) 

(similarly for W B (f) and W' L (j')). 
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